DOC  FILE  COPY  yiDA060827 


NPS55-77-34  ^ 

“WL  FDSTGRADUATE  SCHOOL 


Monterey,  California 


A 


■A  'I 


ii  e* 


S.-5> 


D D C 


NOV  6 1978 


^ I ' ^ PROCEDURE  FOR  ESTIMATING  AN 
jOBJECT'S  ^OSITION  £ASED  ON  TWO 
OR  MORE  BEARINGS  WITH  A ^ROGRAM 
FOR  A _^-59  CALCULATOR  0 '■ 

R.  Neagle/ Forrest  ) (^ / 
September  1977 
( Revi^re^(^/^gi,ia»  iti  7 s >7 


in; 

A 


^ ; 


i 


Approved  for  public  release;  distribution  unlimited. 


o y 


J 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


A 


SeCUWtTY  CLASSIFICATION  OF  THIS  PACE  (Whmn  D»tm  Bnffd) 

REPORT  DOCUMENTATION  PAGE  befoIe'^com^le^Sg^form 

<.  REPORT  number  |2.  GOVT  ACCESSION  NO.  »■  RtCIPIENT'S  CATALOG  NUMBER 

NPS55-77-34^ 


t.  REPORT  number 


4.  title  (mnd  Sublllle)  S.  TYPE  OF  REPORT  & PERIOD  COVERED 

A Procedure  for  Estimating  an  Object's  Position  Technical 
Based  on  Two  or  More  Bearings  with  a Program 

for  a TI-59  Calculator  c.  performing  org.  report  number 


T.  AUTHORf*; 

R.  Neagle  Forrest 


• . CONTRACT  OR  GRANT  NUMBERfa) 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


10.  PROGRAM  element,  PROJECT,  TASK 
AREA  « WORK  UNIT  NUMBERS 


Naval  Postgraduate  School 
Monterey,  California  93940 


|ll.  CONTROLLING  OFFICE  NAME  and  ADDRESS  REPORT  DATE  September  1977 

Naval  Postgraduate  School  (Revised  August  1978) 

Monterey,  California  93940  is.  number  of  pages 

42 

U.  monitoring  agency  name  & ADDRESSCIf  rfl/Zaranl  /ram  ControHIng  0»le»)  IS.  SECURITY  CLASS,  (of  (hla  raportj 

Unclassified 

iSa.  declassification/ downgrading 

SCHEDULE 


is.  distribution  statement  (of  (Ait  Ropetl) 

Approved  for  public  release;  distribution  unlimited. 


I IT.  DISTRIBUTION  STATEMENT  (of  Iho  aSairaci  onfotod  tn  Block  30,  If  dUforonl  from  Ropotf) 


flS.  SUPPLEMENTARY  NOTES 


J? ^ 


1 19.  KCY  words  (ConCInu*  on  rovorco  «(cfo  If  nocooooTY  ond  fdonllfy  by  block  numbor) 


/£>  'Krv  ^ SP 


Navigation  Probability  Area 

Navigational  Error  Elliptical  Probability  Region 

Position  Estimation  Elliptical  Confidence  Region 

Direction  Finding  TI-59  Calculator 

Bearings  Only  Fixing 

CBSTRACTfCo(irJmia'o«rra»araaa73al7nacaaaanr»id73anrIfrsjrMoc*fitaii!aO 

The  Report  provides  a procedure  for  estimating  an  object's  position  based 
on  bearings  taken  from  or  on  the  object  for  two  or  more  stations.  The 
Report  also  provides  a program  for  the  TI-59  calculator  to  Implement 
the  procedure. 


DO  1473  EDITION  OF  I NOV  SS  IS  OBSOLETE 


S/N  0102'014- 6601  I 


security  CLAtSIFICATlOM  OF  THIS  PAGE  Dpi*  Bnlprpdj 


In  this  revision,  the  TI-59  program  has  been 
modified  so  that  a user  may  now  revise  a location  estimate 
by  entering  additional  bearing  data. 

An  example  of  the  use  of  this  option  is  given 
on  Page  16. 


TABLE  OF  CONTENTS 


Page 

I Introduction  1 

II  User  Instructions 

III  Program  Listing 20 

IV  A Development  for  the  Procedure 25 

References 


The  programs  in  this  report  are  for  use 
within  the  Department  of  the  Navy,  and  they  are 
presented  without  representation  or  warranty  of 
any  kind. 


f 


A PROCEDURE  FOR  ESTIMATING  AN  OBJECT'S  POSITION 


! 


BASED  ON  TWO  OR  MORE  BEARINGS  WITH  A 
PROGRAM  FOR  A TI-59  CALCULATOR 

I . Introduction 

A procedure  for  estimating  an  object's  position  with 
bearings  taken  on  or  from  two  or  more  stations  is  developed 
in  Section  IV  of  this  report.  In  the  development  of  the  pro- 
cedure, the  following  things  are  assumed;  The  object  and  the 
stations  are  fixed  on  the  surface  of  a flat  earth  and  the 
position  of  each  station  is  known.  The  error  in  the  bearing 
taken  on  or  from  a station  is  a normal  random  variable  with 
a known  standard  deviation  e and  a mean  of  zero  (if  bias 
exists,  it  is  known  and  removed);  and  station  bearing  errors 
are  independent.  The  user  instructions  for  a TI-59  program 
to  implement  the  procedure  are  given  in  Section  II,  and  the 
program  listing  is  given  in  Section  III. 

As  an  example  to  illustrate  a use  of  the  program,  suppose 
bearings  are  taken  on  an  object  from  three  stations  (1,  2 and  3) 
as  illustrated  in  Figure  1.  Also,  suppose  that  the  assumptions 
stated  above  are  satisfied  and  that  an  initial  estimate  of  the 
object's  position  is  made  and  that  it  is  relatively  near  the 
object.  This  assumption  is  discussed  in  Section  IV. 
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= 10,000  meters. 


r2  = 15,000  meters, 


r^  = 12,000  meters. 


= 38' 


^2  = 3460 
= 127“  . 


Use  of  the  position  estimation  program  with  this  data 
gives  a final  position  estimate  (fix)  determined  by; 

X = -512  meters 
y = - 75  meters 

where  x is  its  East-West  distance  and  y is  its  North-South 
distance  from  the  initial  position  estimate.  The  East-West, 
North-South  xy-coordinate  system  with  its  origin  at  the  initial 
estimate  is  shown  in  Figure  2.  So  the  final  position  estimate 
is  512  meters  to  the  West  and  75  meters  to  the  South  of  the 


initial  position  estimate. 


North 


Initial  Estimate 


• Final  Estimate 


FIGURE  2.  The  Location  of  the  Final  Position  Estimate  with 
Respect  to  the  Initial  Position  Estimate. 


Minimum  area  elliptical  confidence  regions  for  an  object's 
position  can  also  be  found  by  using  the  TI-59  progreun.  The 
centers  of  the  regions  are  at  the  fix,  and  their  axes  lie  along 
the  X*  and  y'  axes  of  the  coordinate  system  obtained  by 
rotating  the  East-West,  North-South  xy-coordinate  system  with 
origin  at  the  fix  through  an  angle  Y . The  angle  Y is  defined 
so  that  it  is  positive  for  a rotation  in  the  counterclockwise 
direction. 

With  the  data  from  the  above  example,  the  program  gives 
Y = -31“ ; so,  the  x'  axis  is  directed  31“  South  of  East.  For 
a confidence  region  with  minimum  area  and  a confidence  level  of 
.9000,  the  ellipse  bounding  the  region  has  a semi -major  axis 
of  2064  meters,  and  a semi -minor  axis  of  1453  meters.  The  area 
of  the  region  is  9.43  square  killometers  or  2.75  square  nautical 
miles.  The  region  is  shown  in  Figure  3. 


North 


FIGURE  3.  A .9000  Confidence  Region  for  an  Object's  Position. 
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In  the  example  discussed  above,  the  position  of  the 
initial  estimate  is  an  input  to  the  program.  If  this  is 
not  desirable,  the  program  can  be  used  to  determine  a po- 
-^'ition  for  the  initial  estimate.  The  position  is  the  in- 
tersection of  the  two  bearing  lines  corresponding  to  the 
first  two  bearings  entered  in  the  program.  Both  options 
are  illustrated  in  Section  II. 

Since,  in  general,  the  smaller  the  bearing  errors, 
the  more  likely  that  the  initial  estimate  will  be  relatively 
near  the  object;  small  bearing  errors  can  be  considered  to 
be  a condition  on  the  use  of  the  procedure. 

Note,  if  the  length  of  the  base  line  joining  the 
first  two  stations  is  small  enough  and  their  bearing  errors 
are  large  enough,  observed  bearing  lines  from  the  two 
stations  may  not  intersect.  If  they  do  not  intersect,  the 
initial  estimate  determined  by  the  program  will  be  at  the 
intersection  of  the  reciprocal  bearing  lines,  and  a gross 
error  can  result. 
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II.  User  Instructions 


The  TI-59  program  to  which  the  user  instructions  in  this 
section  apply  can  be  used  to  calculate  the  quantities  described 
in  Section  I. 

The  program  requires  the  following  inputs: 

1.  the  observed  bearing  from  or  on  an  object  for  two  or  more 
stations ; 

2.  station  positions  relative  to  a reference  position;  and 

3.  the  bearing  error  (standard  deviation)  for  each  observed 
bearing . 

Station  positions  can  be  specified  in  either  of  two  ways. 

In  the  first  way.  Mode  A,  each  station's  position  is  specified 
in  terms  of  its  bearing  a and  its  range  p from  a reference 
position.  In  the  second  way.  Mode  B,  each  station's  position 
is  specified  in  terms  of  its  East-West  distance  x (plus  for 
East)  and  its  North-South  distance  (plus  for  North)  from  a 
reference  position.  The  reference  position  can  be  any  convenient 
location.  For  example,  if  it  were  at  a station,  then  for 
that  station  a = 0 and  p = 0 or  x = 0 and  y = 0. 

The  program  also  requires  an  initial  estimate  of  the  object's 
position.  The  user  has  two  options: 

1.  Let  the  program  provide  an  estimate,  or 

2.  Provide  one  with  the  input  data. 

For  Option  l,the  initial  estimate  is  at  the  intersection  of  the 
bearing  lines  determined  by  the  first  two  observed  bearings  entered 
into  the  program.  For  this  reason,  if  this  option  is  chosen,  the 
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first  and  second  groups  of  data  entered  should  correspond  to 
the  two  stations  estimated  to  have  the  smallest  products 
Although  in  this  option  the  reference  position  cannot  be  at 
the  initial  estimate,  it  can  be  at  one  of  the  stations.  If 
only  two  stations  are  involved,  the  final  estimate  is  at  the 
intersections  of  the  bearing  lines.  (If  the  second  option  of 
either  mode  is  used  with  an  initial  estimate  which  is  not  at 
the  intersection  of  the  two  bearing  lines,  the  coordinates  of 
the  final  estimate  will  differ  from  coordinates  of  the  inter- 
section to  the  degree  of  the  approximations  involved  in  the 
estimation  procedure . ) 

Two  ways  of  providing  confidence  (probability)  region  data 

are  available.  In  the  first  way.  Mode  C,  the  confidence 

(probability)  p is  specified.  In  the  second  way.  Mode  D,  the 

multiplier  k is  specified  where  ko^,  and  ko^',  are  the 

X y 

semi-axes  of  the  bounding  ellipse. 

The  values  of  various  quantities  calculated  by  the  program 
are  either  stored  in  registers  or  appear  in  the  display.  If  a 
PC-IOOA  printer  is  used,  some  of  these  values  will  be  printed. 

The  location  of  calculated  values  and  the  printing  format  is 
given  after  the  user  instructions.  Those  quantities  which  are 
not  described  in  Section  I are  described  below  in  the  User  In- 
structions or  in  Section  IV. 

All  angles  required  or  calculated  by  the  program  are 
in  decimal  degrees. 
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i step 


Instructions 


1. 


If  the  calculator  has  been  in  use  and 
flags  have  been  set  or  the  memory 
repartitioned,  turn  the  calculator  off 
eind  then  on. 


2 . Read  Side  1 and  Side  2 of  Card  1 . 

3.  Read  Side  3 of  Card  2. 


Enter  Press 


MODE  A:  Station  Locations  Specified  in 
Terms  of  Bearing  auid  Range  from 
a Reference  Point. 


4a.  If  the  initial  position  estimate  will 
be  determined  by  the  program,  go  to 
Step  7a.  See  the  note  on  Page  10. 


5a. 

Enter 

the 

initial  estimate’s  bearing. 

a* 

A' 

6a. 

Enter 

the 

initial  estimate's  range. 

P* 

R/S 

7a. 

Enter 

the 

measured  bearing  on  the 

e. 

1 

A 

object  from  a station  or  the  reciprocal 
of  the  measured  bearing  on  a station 
from  the  object. 


8a. 

Enter 

the 

station's  bearing. 

“i 

R/S 

9a. 

Enter 

the 

station's  range. 

Pi 

R/S 

10a. 

Enter 

the 

bearing  error. 

®i 

R/S 

11a.  Repeat  Steps  7a,  8a,  9a  and  10a  for  all 
stations.  The  number  of  repetitions  i 
appears  in  the  display  after  Step  10a. 


8 


1 


Display 


i 


1 


Step 

Instructions 

Enter 

Press 

Display 

MODE 

B:  Station  Locations  Specified  in  Terms 
of  East-West  Distance  and  North-South 
Distance  from  a Reference  Point. 

4b. 

If  the  initial  position  estimate  will 
be  determined  by  the  program,  go  to 

Step  7b.  See  the  note  on  Page  10. 

5b. 

Enter  the  initial  estimate's  East-West 
distance . 

X* 

B' 

6b. 

Enter  the  initial  estimate's  North- 
South  distcince. 

y* 

R/S 

7b. 

Enter  the  measured  bearing  on  the 
object  from  a station  or  the 
reciprocal  of  the  measured  bearing 
on  a station  from  the  object. 

e. 

1 

B 

8b. 

Enter  the  station's  East-West 
distance. 

R/S 

9b. 

Enter  the  station's  North-South 
distance. 

R/S 

10b. 

Enter  the  bearing  error. 

®i 

R/S 

i 

11b. 

Repeat  Steps  7b,  8b,  9b  and  10b  for  all 
stations.  The  number  of  repetitions  i 
appears  in  the  display  after  Step  10b. 

BOTH 

MODES 

12. 

Calculate  the  East-West  distance. 

R/S 

the  North-South  distance,  the 

bearing  and  the  range  of  the 

position  estimate  relative  to  the 

reference  position.  Also  calculate 

the  rotation  angle  y , and  the 

standard  deviations  Ca.  and  a^,. 

x'  y' 

To  include  additional  bearing  measure- 
ments after  this  calculation,  gc  to 
Step  18. 
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13.  For  confidence  (probability)  region 
calculations,  go  to  Step  14  if  the 
confidence  (probability)  for  the 
region  is  specified.  If  )c  is 
specified  where  ko-.  and  ko.. 

X y ' 

are  the  semi-axes  of  the  bounding  ellipse 
with  the  larger  the  major  axis,  go 
to  Step  16. 

14.  Enter  p,  the  confidence  level  p C Area 

(probability)  and  calculate  k, 

ka»,,  and  the  area  of  the  region. 

X y 

(The  area  units  correspond  to  the 
distance  units  used.) 

15.  For  a different  value  of  p,  go  to 
Step  14. 

16.  Enter  k and  calculate  the  confidence  k D Area 

level  (probability)  p,  ko..,  ko,.  and 

X y 

the  area  of  the  region.  (The  area 
units  correspond  to  the  distance  units 
used . ) 

17.  For  a different  value  of  k,  go  to 
Step  16. 

18.  To  include  an  additional  bearing  measure- 
ment from  either  a new  or  old  station, 

go  to  Step  7a  if  using  Mode  A or 
Step  7b  if  using  Mode  B. 


NOTE:  If  a data  entry  error  occurs  in  either  mode, 

press  RST  and  then  use  the  following  pro- 
cedure: For  Option  1,  return  to  Step  7 
and  repeat  all  data  entries.  For  Option  2, 
return  to  Step  5 and  repeat  all  data  entries. 

Also,  if  a position  estimate  is  to  be  de- 
termined for  a new  object  position  or  if 
a new  mode  is  to  be  used,  follow  this 
instruction . 
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NOTES ; 

a)  The  program  printing  format  is  given  below: 

For  the  initial  data,  i = 1,  2,  ...  , n with  one  space 
between  groups: 


Mode  A 


a* 


«i 


initial  estimate 
if  provided 


Mode  B 

X* 

y* 


e. 

1 


e . 
1 


The  format  for  the  calculated  position  data  is: 

X 

y 


a 

P 


For  the  confidence  (probability)  region  portion  of  the  program 
the  format  is  after  pressing  either  C or  D; 

P 

k 


kdA , 

semi -axis 

ko^  , 

y' 

semi -axis 

Area 

b)  The  following  data  is  stored  in  the  indicated  registers; 


Data 

Registers 

X* 

R38 

y* 

R39 

Y 

R29 

X 

R30 

y 

R31 

a 

R32 

P 

R33 

a 

x' 

R16 

0 ^ , 
y 

R17 

P 

R14 

k 

R15 

ka-'  , 

R18 

kPr;, 

R19 

y 

12 

I 

J 


Four  data  tapes  for  a sample  problem  are  given  below. 


Distance  units  have  not  been  specified,  but  they  could  be 
meters  for  example.  Angles  are  in  degrees.  Option  1 (initial 
estimate  not  provided)  for  Mode  A and  Mode  B is  indicated  by 
A and  by  B and  Option  2 (initial  estimate  provided)  is  indicated 
by  A ' and  B ' . 

For  each  mode  and  each  option,  the  input  data  are 
indicated.  The  data  determine  the  relative  locations  of  three 
stations  as  well  as  the  observed  bearing  of  an  object  from 
each  station. 

For  A and  B,  the  reference  location  is  at  Station  1 
and  the  initial  position  estimate  (determined  by  the  program) 
is  at  the  intersection  of  the  bearing  lines  for  Station  1 and 
Station  2. 

The  intersection  has  coordinates  x*  = 906.4853528 
and  y*  = 17296.77092  with  respect  to  Station  1. 

For  A'  and  B',  both  the  initial  estimate  and  the 
reference  location  are  at  the  intersection  of  the  bearing  lines, 
so  a*  = 0 and  p*  = 0 and  x*  = 0 and  y*  = 0. 

The  data  for  A,  B,  A'  and  B'  are  all  equivalent, 
and  each  solution  gives  the  same  data  for  a confidence  (prob- 
ability) region  calculation.  A tape  with  confidence  (prob- 
ability) region  results  for  both  Mode  C and  Mode  D which 
correspond  to  A,  B,  A'  and  B'  is  given  with  the  first  four 
data  tapes. 
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£.  ('  Zf  m 

02 

4. 

1 0000. 

P2 

3 m 

®2 

O • J»  a 

213. 

303. 

03 

20000. 

“3 

3. 

1 4000. 

P3 

*m.*  a 

®3 

303. 

129. 

5867755 

573. 5078933 

X 

8717 

. 797886 

1 b4t  ,3.  7 i 233! 

y 

8. 

1. 995471725 

o 

18472. 70157 

p 

8974567 

-834. 

058683:5 

-7. 325492245 

y 

201 . 

75840 1 9 

787,  3663755 

8 9 8 . 

03931  1 1 

1233. 080777 

Ca 

y 

5392259 

-7r. 

1 1*  a 

3663757 

1233 

. 080776 

THIS  PAGE  IS  BEST  QUALITY  PRACIlCABlH 
raOM  COPY  f iO  DUG  - 


B 


B' 


C or  D 


3. 
0. 
0. 

4. 


33.  02 

■9986.  295348  ^2 
523.  3595624  72 
3.  ©2 


e. 


303. 
7624. 94649 
11741.38795  73 

8. 


573 


3 

«3 


:iS7«y37  X 
16462.71223  7 
1.995471723  ® 
16472.70157  p 

-7.  325392245  y 
787.  3663756 
1233.080777 

7* 


U. 

0. 


■906.  4S5352S 
■l729t..  77092 

4. 


-10892. 7807 
■16773.  41136 


303. 

67 1 8.  46 1 1 37 
-5555.  382969 


-332. 8974589 
-834. 0586906 
201. 7584019 
898. 0393184 

-7. 325392239 
787. 3663755 
1 233.  U 8 0 1'  i'  1' 


X* 

Q.  9 

P 

V* 

2.  1 4 5 '7‘t^  01.126 

k 

Jr 

1639. to  1493 

9l 

2646. 149454 

ka^, 

14046  664.  97 

Ar)^a 

. 3 64664^1 6 8 

P 

k 

®2 

1574. 732751 

ka^. 

2486. 161553 

ko^. 

X2 

^2 

«2 

12200517. 6 

63 

X3 

^3 

63 


7 

a 


Y 

X' 

7' 


irilS  ?AG£  IS  BEST  QUALITY  PT^AP.TTrJLRT^E 
EROll  COrY  rUSiilSiiiii)  lODDO 


The  following  data  tape  illustrates  the  effects  of 
using  only  the  bearings  for  the  first  two  stations.  The  tape 
is  for  Mode  A,  Option  1.  The  values  stored  in  Registers  38 
and  39  (the  x and  y coordinates  of  the  intersection  of  the 
bearing  lines  from  Station  1 and  Station  2 with  reference, 
in  this  case,  to  Station  1)  are  also  listed  on  the  tape  (as 
well  as  given  above).  And,  as  can  be  seen,  the  initial 
estimate  aind  final  estimate  correspond. 

The  data  tape  also  illustrates  the  use  of  additional 
bearing  data  to  revise  a position  estimate.  The  data  for 
Station  3 printed  after  the  first  confidence  (probability) 
region  calculation  results  was  entered  by  again  repeating 
Step  11a,  and  the  remaining  results  were  obtained  by  next 
repeating  Step  12  and  then  Steps  14  and  16.  Note,  these 
results  are  the  seune  as  the  corresponding  results  for  Mode  A 
on  page  14. 
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To  obtain  the  results  given  in  Section  I,  use  A* 
and  take  the  reference  position  at  the  initial  estimate 
(a*  = 0,  p*  = 0).  Then  = 218°,  a2  = 166°  and 
^‘3  307°.  The  data  tape  for  the  calculation  is  given  below. 
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18 


A' 

0. 

a* 

0. 

P* 

35. 

2 IS. 

“1 

1 uuuu. 

Pi 

4. 

351. 

02 

1 66. 

02 

1 5000. 

P2 

1 • 

131. 

03 

307. 

“3 

12000. 

P3 

5. 

®3 

-511. 961856 

X 

-75.  43753883 

y 

261. 617789 

a 

517. 4898687 

P 

■31. 23492683 

Y 

677. 2632305 

0.  , 

961. 6888632 

0^ 

y* 

0.  9 

p 

2.  145966026 

k 

1453. 383883 

ko^ 

2063. 751628 

ka*' 

9422966. 381 

Ar^; 

. 86-^6647168 

P 

k 

1 354. 52646 1 

ko  ^ 

1 Q O • ^ "T  •* 

1 -•  B O l'  1 i t' 

ko*; 

8184684. 605 

y ' 

Area 

1 

Ill . Program  Listing 

Before  entering  the  program,  press  2nd  and  then  CP  or 
turn  the  calculator  off  and  then  on.  Next  enter  5 in  the  display, 
press  2nd  and  then  Op  17.  This  repartitions  the  calculator's 
memory  so  that  the  complete  program  can  be  entered. 

Before  recording  the  program,  enter  6 in  the  display, 
press  2nd  and  then  Op  17.  This  returns  the  calculator's  memory 
to  the  normal  partition  (479.59).  Returning  the  calculator  to 
the  normal  partition  allows  the  two  program  cards  to  be  read 
in  the  normal  partition  without  forcing.  When  the  program 
is  used,  it  repartitions  the  calculator  so  that  Bank  3 
registers  are  program  registers. 
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42 
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19 
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42 
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19 

D* 

173 

O'? 

1 
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37 
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75 

- 
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94 
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oo 

w w 
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22 
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87 
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1 89 
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01 
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01 

01 

191 

93 

93 
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n* 
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75 

- 
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43 
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39 

39 
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95 
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22 
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205 
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0 1 0 1 
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216 
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222 
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^ 
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259 

44  SUM 

260 

21  21 
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IV.  A Development  for  the  Procedure 


I 


1 


In  the  development  for  the  estimation  procedure  given  here, 
all  angles  are  in  radians  and  the  assumptions  stated  in  Section  I 
apply. 

Figure  4 shows  three  bearing  lines  from  the  ith  of  n stations. 


One  IS  the  observed  bearing  line 


of  an  object.  One  of  length  goes 

to  the  origin  of  an  xy-coordinate 
system  located  at  the  object's  un- 
known position.  And  one  of  length 
goes  to  an  initial  estimate 
with  known  position  but  unknown 
coordinates  (x,y) . Note,  esti- 
mates for  -X  and  -y  estimate 
the  object's  position.  To  find 
estimates  -x  and  -y,  consider 
the  arc  coordinates  u^  = ^^(6^-4)^) 
of  the  observed  bearing  line  and 
Vj^  = of  the  bearing  line 

to  the  point  (x,y) . They  are  de- 
fined by  the  three  bearing  lines 
and  the  circle  of  radius  which 

goes  through  the  object's  position 
and  which  is  centered  on  the  station 
as  shown  in  Figure  4. 


FIGURE  4.  Problem  Geometry. 


25 


t 


By  defining  (all  angles  in  radians), 

u.  = V.  + w..  Note,  0.-S.  is  known,  but  is  not. 

However,  can  be  expressed  in  terms  of  x and  y,  and,  to  first 

order,  = x cos  6^  - y sin  0^;  so,  if  tan  =*  (3^-(j)^)  for 

i = l,2,...,n,  that  is,  if  (x,y)  is  relatively  near  the  object's 
position,  + x cos  6^  - y sin  6^  for  i = l,2,  — ,n. 

In  this  development,  the  initial  estimate  is  assumed  to  be 

relatively  close  enough  to  the  object's  position  so  that  the  above 

approximation  for  u^  can  be  used  and  so  that  can  be  taken 

equal  to  r^.  With  this  assumption,  u^  is  a function  of:  6^, 

the  observed  value  of  a random  quantity;  the  known  parameters  r^ 

and  and  the  unknown  parcuneters  x and  y. 

If  a distribution  for  the  0^  can  be  specified,  then  a 

distribution  for  the  U.  can  be  determined  and  standard  estimates 

1 

A 

X and  y for  x and  y can  be  considered.  Maximum  likelihood 

estimates  are  discussed  in  this  section.  Each  0^  is  taken  to 

2 

be  a normal  random  varieJale  with  mean  and  variance  e^.  And 

the  n random  variables  0^,  i = l,2,...,n  (one  for  each  station) 

are  taken  to  be  independent. 

The  likelihood  for  a sample  6,,  6-,  ...  , 0 is  then 

12  n 


V - ; (;i^)  [-  5 f ‘“rh''/®!] 

and  the  likelihood  for  a corresponding  sample  Uj^ , u^  is 
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,(u,  ,u„,  . . . ,u  ) = n ( — \ exp  f-  y z 

^ ^ " i=l  \ /2F  a.  / L ^ 1 1 


where  a.  = r.e.  (with  e.  in  radians)  since  u.  = r.(0.-i|).). 

Ill  1 1111 

By  definition,  the  maximum  likelihood  estimates  of  x 

/S  A 

and  y are  the  estimates  x and  y which  make  L (Uj^,U2,  . . . ,u^) 

a maximum.  In  this  case,  making  L (u^^  ,U2 , . . • , u^)  a maximum  is 

" n 2 2 ^ 

equivalent  to  making  a minimum.  So,  to  find  x 

A 

and  y,  solve  the  following  two  equations  for  x and  y: 


9()ln  L) 


= 0 and 


9(«,n  L) 


= 0 . 


The  solutions  are  x = x and  y = y/  x and  y are  the 

maximum  likelihood  estimates.  With  w.  = r.(0.-3.)  and 

1 111 

the  conditions  assumed  above  these  two  equations  are  linear 
equations  in  x and  y.  And, 


r 2 

I [w.  + X cos  3-  - y sin  3-]  (cos  3 . ) /o . = 0 

1 ^ X XX 

? 2 

I [w.  + X cos  3.  - y sin  3 ] (sin  = 0 

^ X X X XX 


And,  in  terms  of  the  following  quantities: 

2 2 7 

A = 1 (cos  3^)/c^  , B = Z(sin  3^^  cos  3^)/a^  , 


C = E(sin^  3^)/o?  , 


E = E(w^  sin  3^)/af  , 


D = E (w^  cos  , 


the  equations  are : 


Ax  - By  = -D 

Bx  - Cy  = -E  . 

So  the  solutions  are: 

X = (BE  - CD) /(AC  - B^) 

y = (AE  - BD)/(AC  - B^)  . 

A confidence  region  can  be  constructed  about  an  estimated 
position.  In  order  to  indicate  how  this  is  done,  a probability 
region  about  the  true  position  will  be  considered  first. 

Note,  X and  y are  values  of  random  variables.  If  a 
new  set  of  bearings  0 , 62 » . • • » 6^^  is  observed  (for  a fixed  initial 
estimate  and  object) , in  general,  a new  pair  of  values  x and 
will  be  obtained. 

/V  A. 

If  X and  Y represent  these  random  variables,  then 

/V  1 ^ 

X = 5—  y (VI. /a.)  (B  sin  0.  - C cos  0.) 

(AC-B*^)  1 ^ ^ ^ 


Y = =— *—  y (W./0-)  (A  sin  0.  - B cos  0.) 

(AC-B^)  1 ^ ^ ^ 

with  = r^(0^-0j^).  (W^  is  the  random  distance  intercepted 

along  the  ith  arc  between  the  bearing  lines  defined  by  0^  and 

6..) 
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4 


►<> 


Note,  X and  Y have  a bivariate  normal  distribution. 


since  they  are  a linear  combination  of  the  n normal  random 
variables  ,W2 , . . . ,W^,  or  equivalently  of  the  n normal 
random  variables  Q.,©,,...,®  . Also  E (W . ) = r.  (<j).-8.). 

If  6^  = for  i = l,2,...,n,  that  is,  if  the  initial 
estimate  of  the  object's  position  is  at  the  object's  position, 

E(Wj^)  = 0 for  i = 1,2, ...,n.  And,  therefore,  E(X)  = 0 and 
E(Y)  =0.  So,  in  this  case,  the  "location"  of  the  bivariate 

A A 

normal  distribution  of  a point  (x,  Y) , the  random 
coordinates  of  the  object's  estimated  position,  is  the  same  as 
that  for  the  point  (-X,  -Y)  and  both  are  centered  on  the 

object's  position.  However,  the  "location"  of  the  distribution 

/\  /v 

of  (-X,-Y)  is  independent  of  the  location  of  the  initial 
estimate  when  the  coordinates  (-X,-Y)  refer  to  a coordinate 
system  with  origin  at  the  initial  estimate.  This  fact  simplifies 
the  establishment  of  a confidence  region  about  the  location  of 
an  estimated  position. 

A region  of 'minimum  area  for  a given  probability  of  con- 
tainment of  an  estimated  position  can  be  determined.  The  region 
is  bounded  by  an  ellipse  which  is  centered  on  the  object's  position 
and  whose  axes  lie  along  the  axes  of  an  x 'y ' -coordinate  system 
obtained  by  rotating  the  xy-coordinate  system  centered  on  the 

object's  position  through  an  angle  y . In  this  system,  a = 0, 

^ xy 

that  is,  X'  and  Y'  are  independent  normal  random  variables. 
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The  two  coordinate  systems  are  illustrated  in  Figure  5. 
The  coordinates  of  a point  in  the  two  systems  are  related  by 
y'  y x'  = X cos  Y + y sin  y 

y*  = -X  sin  Y + y cos  y 

X 

These  relations,  along  with 

^ =0,  imply: 

x'y' 

FIGURE  5.  Rotation  Geometry. 


where  y,  the  angle  of  rotation  of  the  coordinate  axes,  is 
positive  in  the  counterclockwise  direction. 

With  the  initial  estimate  of  the  object's  position  at 
the  object's  position  (3^^  = i = l/2,...,n),  so  E(W^)  = 0 
and  Var(W^)  = a?, 


I 


I 


[ 

f 

I 

f 


■■  ~ ^ 

J 

! 

] 

- i-5-5-  I a/ah  (B  sin  6.  - C cos  3.)^,  ! 


- « I (1/a?)  (A  sin  3.  - B cos  3^)^ 
y (AC-B"^)  1 ^ ^ ^ 

and 

1 ^ 2 

O/s/s  = 5— *■  'l  a/a. ) (B  sin  3 . - C cos  3 . ) (A  sin  3 • - B <30%  3 ■ ? . 

xy  (AC-B^)"^  1 ^ ^ ^ ^ ^ 


Using  the  definition  for  A,  B and  C,  the  above  become' 


7~> 

X (Ac-e  ) 

y (AC-B'^) 

and 

r,  - B 

^ ~ 5~  • 

xy  (AC-B^) 

So,  tan  2y  = 2B/(C-A)  for  3^  = <|)^,  i = l,2,...,n. 

With  the  object's  position  known  and,  hence,  (j)^  known 

2 2 

for  1 = l,2,...,n,  the  above  equations  for  a , a , a and  y 

^ /S  A 

X y xy 

can  be  used,  since  the  initial  estimate  of  the  object's  position 

can  be  teJcen  as  the  object's  position. 

With  values  for  a^,  a^,  and  Y,  values  for 

X y xy  x' 

and  can  be  found  by  using  the  equations  in  the  middle  of 

y' 

Page  30.  And  then,  the  prob^U^ility  that  an  estimated 

position  will  be  within  an  ellipse  of  semiaxes  ko^  and  ka^ 

x'  y ' 
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which  is  centered  on  the  object's  position  can  be  found.  It  is 
i - exp(-)c  /2).  (This  result  follows  from  integrating  the  bi- 
variate normal  density  over  the  ellipse.)  And  the  area  of  the 

ellipse  is  irk^CA  O/s 
X ' y ' 

I Given  estimates  x and  y found  by  using  the  relations 

on  Page  28,  an  ellipse  with  semi-axes  ko^  and  ko^  centered 
! x'  y ' 

on  the  point  with  coordinates  (-x,-y)  in  a coordinate  system 

' with  origin  at  the  initial  estimate  and  oriented  as  indicated 

i 2 

by  Y is  a 1 - exp(-k  /2)  confidence  region.  This  follows 

I from  the  bivariate  normal  distribution  of  -X  and  -Y  which 

in  this  system  is  centered  on  the  object's  position.  The  ellipse 
2 2 

is  defined  if  and  are  known  (the  covariance  matrix 

X y xy 

is  known) . And  to  the  degree  of  the  approximations  involved, 
this  can  be  assumed  to  be  the  case.  In  particular,  by  assuming 
! the  initial  estimate  of  the  object's  position  is  at  the  object's 

I 

1 

position,  which  is  consistent  with  assuming  (3^-<()j^)  is  small, 

2 2 

' values  for  of  , and  Y can  be  obtained  by  using 

X y xy 

the  relations  on  Page  31.  These  values  can  then  be  used  to 
I 2 2 

\ determine  and  by  using  the  relations  on  Page  30. 

1 X ' y ' 

And,  then,  with  a value  for  k,  a confidence  region  can  be 
constructed.  To  the  degree  of  the  approximations  involved,  the 
j shape  of  the  confidence  region  is  independent  of  both  the  object's 

I position  and  of  the  initial  estimate  of  the  object's  position. 


1 
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For  the  case  where  bearings  are  taken  from  the  object 
on  two  or  more  stations,  is  the  reciprocal  of  the 

bearing  taken  from  the  object. 

A discussion  for  this  and  for  other  bearings  only 
position  estimation  procedures  for  situations  similar  to  the 
one  considered  here  is  given  in  Reference  1 listed  below. 
Reference  2 gives  an  equivalent  bearings  only  procedure.  It 
also  gives  a range  only  procedure,  a range  and  bearing  procedure 
and  HP-9830A  progreuns  with  which  to  implement  the  procedures. 
Using  the  fix  determined  by  two  lines  of  bearing  as  the  initial 
estimate  was  suggested  by  this  reference. 


The  equations  used  in  the  program  to  determine  (x*,  y*) , 
the  coordinates  of  the  fix,  are: 


X*  sin 

(62-61) 

= [Pi 

sin 

^“l“®l^ 1 

sin 

®2 

- fP2 

sin 

^“2~®2^  ^ 

sin 

®1 

y*  sin 

(62-61) 

= [Pl 

sin 

^“l”®l^ ] 

cos 

®2 

- [P2 

sin 

(02-62)] 

cos 

®1 
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